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Abstract 

Mixing effect in a stratified fluid is considered and examined. Euler equations for incom- 
pressible fluid stratified by a gravity field are applied to state a mathematical problem and 
describe the effect. It is found out that a system of Euler equations is not enough for a for- 
mulation of correct generalized problem. Some complementary relations are suggested and 
justified. A numerical method is developed and applied for study of processes of vortex de- 
struction and mixing progress in a stratified fluid. The dependence of vortex destruction on a 
stratification scale is investigated numerically and it is shown that the effect increases with the 
stratification scale. It is observed that the effect of vortex destruction is absent when the fluid 
density is constant. Some simple mathematical explanation of the effect is suggested. 



1 Introduction 



The system of Euler equations for an incompressible fluid stratified by a gravity field is consid- 
ered. It is known that internal gravity waves can propagate within fluids stratified by density in a 
gravity field. Internal wave breaking is one of interesting nonlinear effects, which is characterized 
by disintegration of waves with formation of typical spots with intensive small-scale convection 
inside them. The spots are often termed convective or turbulent ones because the convection in- 
side them looks like turbulence. Static stability in the fluid becomes recovered with the course of 
time, but the density inside the spot remains different from the surrounding stratification. The 
phenomenon of internal wave breaking going with formation of a convective spot often is termed 
"internal wave mixing". 

The effect of internal wave breaking is often observed in the ocean. The phenomenon is reg- 
istered in the atmosphere [[5| as well. Broken down internal waves reorganize stratification with 
time. Therefore, it is impossible to understand and explain stratification of the ocean or the atmo- 
sphere without taking into account the effect. Intensive intermixing of water by internal waves 
creates good nutrient medium for plankton in the ocean. Regions of the ocean, where wave break- 
ing happens, are generally places of heightened biological productivity. It stimulates growing 
interest to them and to the physical effect, being cause of heightened biological productivity. 

Many papers are dedicated to experimental study of the phenomenon J6l [Zl El O [TU]]. In McE- 
wan's experimental papers, the phenomenon has been studied by means of thin laser measure- 
ments. Internal waves were excited in a laboratory tank with sizes of 25 cm x 50 cm x 25cm filled 
with a stratified fluid. Stratification is formed by dependence of water saltiness on height. The 
liquid density in the tank varied with height only by 4%. Therefore, it is possible to term this 
stratification as a weak one. The composite experimental equipment including two lasers has al- 
lowed to watch very small-scale density fluctuations originated owing to vortex destruction in the 
stratified fluid and to discover many details of the effect under study. Observation of small-scale 
structures had been based on dependence of refraction index of light on liquid density. McEwan 
has discovered and studied different stages of evolution of wave breaking and mixing j9|, [II 01 ]: 
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overturning, development of interliving microstructure, restoring of static stability. According 
to experiments |9|, [1 1 1| . the initial phase of evolution of a nonlinear wave finishes formation of 



tongues of more heavy liquid inside more light, and vice versa. Further, the effect of overturn- 
ing develops, and tongues disintegrate into waves of much smaller scales. The static stability is 
restored in course of time, but the fluid stratification becomes changed. 

Let us discuss reasons of formation of tongues. The velocity field within a local vortex is 
inhomogeneous, and top speeds are achieved at some distance from the vortex center. It leads 
to formation of the tongue of a heavy liquid, permeating into strata of a light liquid, and of the 
tongue of a light fluid, permeating into strata of a heavy liquid. Different fluid particles cover 
different path lengths for the same time, and difference of path lengths increases with time. If 
we neglect influence of gravity, such fluid movement within a vortex of an elliptic shape leads 
to distortion of isopycnic lines (lines of equal density) and to reeling of isopycnic lines onto an 
ellipse. The ellipse semiaxes are determined by distances from the vortex center up to particles 
having top speeds in the vortex. Some thin layer structure in which stratums of light and heavy 
fluids are alternating is an outcome of local rotating fluid movement. Imposed gravity field gives 
birth to wave processes, and it entails breaking down of the regular thin layer structure. So, we 
observe mixing of the fluid. 

Due to the phenomenon mechanism, whatever weak the stratification would be, the density 
gradients will become arbitrary large with the course of time. Therefore, when we start math- 
ematical modeling of the phenomenon, it is natural to presume a solution of nonlinear fluid 
equations may be nonsmooth and may contain discontinuities. It seems to be of interest to study 
nondifferentiable, generalized solutions of hydrodynamic equations of a stratified incompressible 
fluid placed into a gravity field. 

Further, the fluid is supposed to be an ideal one. Estimates show that dissipative effects are 
important at late stages of evolution. Therefore, the nondissipative problem is of certain inter- 
est. To be sure, taking into account of dissipative terms should increase solution smoothness. 
The smoothing action of dissipation is explained by inhibition of waves with scales smaller than 
some minimally possible scale dependent on dissipative constants. At numerical modeling, the 



smoothing action of dissipation shows one's worth only if the grid step is less than the minimally 
possible scale. Often it is not advantageous or impossible to increase resolution up to a minimally 
possible scale. It is one more argument for study of a nondissipative problem. 

Processes of mixing in viscous two-layer 3D- fluids with unstable stratification were numeri- 
cally explored in [1111] , [1121] (see also the literature quoted there). Authors for raising productivity 
of computations have utilized parallelizing of calculations. The present paper essentially differs 
from [II 1| . [II 2] because in mentioned papers the instability develops at the expense of unstable 
stratification (a heavy liquid has been initially placed above a light one). In the present paper, 
the stratification is continuous and stable; but in time the instability develops in itself, that is, 



under other mechanism than in [II II] , [112] . The instability is shaped by an initial vortex as it 



occur in occur in experiments of Q, [I10IJ . Furthermore, the problem under consideration is a 



two-dimensional one, and dissipative processes are not taken into account, excepting (3.6). As 
against the mentioned papers, we are interested not only in simulation outcomes, but also in sub- 
stantiation of mathematical statement of the problem and in admissible numerical methods. It 
is explained by nontrivial items in the problem. Researched questions are important for viscous 
liquid as well. 



Monographies HI, and papers HSJ, US], HZ]], dS], [Qg, $20§ are dedicated to mathe- 



matical research of linear equations for small-amplitude internal waves. As for nonlinear Euler 
equations for incompressible fluid, correctness of many problems is not proved yet [121] . Multi- 
mode approach [2], in a finite-mode model directly related to McEwan experiment demonstrates 
important role of nonlinear dispersion [[3| We study nonlinear equations and show that a sys- 
tem of Euler equations is not enough for statement of a correct generalized problem in case of a 
stratified fluid. The effects of wave overturning and fluid mixing are simulated numerically. 

The study of dependence of phenomenon of vortex destructions on a stratification scale is 
another purpose of the research. By means of numerical experiments, we shall demonstrate the 
decreasing of overturning effects with diminution of stratification scale H. On the contrary, the 
effect grows with increasing stratification scale H. Nevertheless, it have been discovered that if 
we take a limiting case when the fluid density is a constant, the effect of vortex destruction is 



absent. 



2 Basic system of equations 

We suppose a fluid behavior is described by Euler equations for an incompressible fluid 

dp dpu dpw g dpu dpu 2 dpuw dp (2 1) 

dt dx dz dt dx dz dx' 

dpw dpuw dpw 2 dp du dw 

dt dx dz dz dx dz 

Here p is the density, p is the pressure, u and w are the horizontal and vertical mass velocities 
of the fluid, t is time, x, z are the horizontal and vertical coordinates respectively, g is the free fall 
acceleration. 

It is supposed as well that the nonperturbed fluid is stratified in density exponentially: 

Po(z) = Pooexp(-^-J, H = 6.23m, p 00 = 1000kg/m 3 . 

The chosen value of H approximately corresponds to the one for McEwan's experiments [[9|, [I10IJ . 
In considered McEwan's experiments, the density stratification was linear, but the stratification 
scale was much greater than the size of the tank; therefore, the linear and exponential stratifi- 
cations are practically indiscernible. An exponential stratification is more convenient for mathe- 
matical research of the equations. 

The boundary condition is (v,n) = at the boundary dCl of the domain Q, where v = (u,w) and 
n is a normal line to the boundary. Keeping in mind McEwan's experiments, we take a rectangle 
with horizontal dimension 50 cm and vertical dimension 25 cm as the domain Q. 



3 Statement of the generalized problem 



3.1 Preliminary analysis 



One of tested methods of numerical integration of Euler equations (2.1 1 was based on Galerkin 



method [1221 ]. The solution was searched for as a series in terms of functions dependent on vertical 
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coordinate z and belonging to a complete orthogonal system of functions, the expansion coeffi- 
cients were dependent on horizontal coordinate x and time t. From Euler equations, by known 
recipe we have deduced the equations for expansion coefficients [3J. We term these equations 
as a set of equations for wave modes. The set of equations for wave modes have been solved 
numerically, and then the velocity and density have been calculated using the series. 

In such a way we successfully simulated an initial phase of evolution of the nonlinear wave, 
which was finished with formation of tongue of a heavy liquid inside a light one, and vice-versa. 
A typical scale of the tongues is about 2cm. Formation of such tongues is (observed) watched 
in laboratory experiments as well. In laboratory experiments, these tongues further disintegrate 
into waves of very small scales. However, we have miscarried this effect, though formally it is 
possible to obtain solution of the problem for all time. 

In-depth study of the problem has given the discovery that it is impossible to simulate cor- 
rectly the nonlinear process under study within the framework of a standard Galerkin method. 
The reason is that partial derivatives do not exist for functions with discontinuities. In such cases 
we should search for the nondifferentiable solution as a generalized, nonclassical one. In 1954, 
Lax, considering quasilinear equations, has shown that nonlinear equations admit definition of 
different, nonequivalent classes of generalized solutions [1231 ], [I24IJ . [1251 ], [I26IJ . (see also [I27|]). 

It is explained by the fact the main equations are crashed at discontinuity surfaces, so we 
need some additional conditions to match solutions on discontinuity surfaces. However, the con- 
ditions surface may be taken different, if we consider the question only from a mathematical 
standpoint. Therefore, convergence and stability of the numerical method yet do not guarantee 
that an obtained generalized solution is valid from the physical point of view. Moreover, only one 
of possible definitions of a generalized solution may fit the physical sense. Sufficient conditions 
of physical validity of the generalized solution are not formulated yet, but necessary require- 
ments are obvious. A generalized solution being adequate to the physical sense, should meet the 
condition of conservation of mass, momentum and energy. 

The considered equations for incompressible fluid differ from surveyed by Lax (in what sense), 
but Lax's argumentation is applicable to them as well. Analysis has shown that the conventional 

6 



Galerkin method does not allow achieving simultaneous conservation of mass, momentum and 
energy for an obtained nondifferentiable solution. May be it is possible to correct somehow con- 
ventional Galerkin method that would allow to simulate correctly solutions with discontinuities. 
However, the authors cannot offer acceptable adjusting of a Galerkin method. At the same time, 
finite-difference methods give major freedom that permits choosing a numerical scheme ensuring 
warranted realization of fundamental conservation laws. 

3.1.1 Energy functionals in the theory of linearized hydrodynamic equations for the strati- 
fied fluid 

First, we remind how the theorems of existence, uniqueness and stability for linear problems 
are proved. Ordinarily linear wave equations possess some quadratic functional conserved over 
solutions. This conserved quadratic functional is called a wave energy functional. For example, 



if we linearize equations (2.1 1 with respect to a background density, we come to known equations 
of the theory of small-amplitude internal waves for which the law of wave energy conservation is 
(particular in [141128]): 

^p = 0, H lin = l - f p (z)[(u 2 + w 2 ) + gH4> 2 ]dn (3.1) 



Here the function <p(x,z,t) is defined by the formula p(x,z,t) = p (z)(l + <p(x,z, t)). In (3.1| the 



term containing velocities means a kinetic energy, and the term containing lgH<p 2 is a potential 



energy. Functional (3.1 1 is of great importance in mathematical theory of equations for small- 
amplitude waves: proofs of theorems of existence, uniqueness, stability for the linearized hydro- 
dynamic equations essentially uses conservation of functional H/,„ over solutions. 
The linearized fluid equations also possess other conservative functional of energy 



« = J Po( z ) \(u 2 + w 2 ) + gz<p 



dCl (3.2) 



We will name this functional a functional of hydrodynamic energy [1291] . In (3.2 1, the term con- 
taining velocities means the kinetic energy and the term containing gz<p is the potential energy. 

Here both conservative energy functionals are written to underline their distinction. Differ- 
ent formulas for energy functional are possible because energy is determined up to a constant. 



If we add any other conservative functional to the energy functional, we obtain a conservative 
functional again, which we may term an energy functional, but the expression for energy density 
becomes other. Therefore, different formulas for energy functional exist. 



We have selected the form of wave energy functional (3.1 1 so that the energy to be zero at lack 



of waves. For this reason, functional (3.1 1 is named wave energy functional. Functional ( 3.2[ ) is 



more known, but not possesses necessary mathematical properties. 
3.1.2 Wave energy functional for nonlinear equations 

Having the purpose to develop a theory for nonlinear equations, it is naturally to extend wave 
energy functional ( 3. 1 [ > to nonlinear case. The following theorem is valid. 



Theorem 1. Let a solution to equations (2.1 ) be differentiate. Then the functional 

,2 , ,„2 



Hnoni = Jjp + 2 W + PS Z+ gH |p ln|-^y J + (p (z) - p) 

D 



dn (3.3) 



is conserved on solutions. The integrand in (3.3) is strictly non-negative at p > 0. The functional (3.3) 



turns into the functional (3.1 ), when we take the small- amplitude limit. 



Proof. Conservation of the functional ( |3.3| follows the fact that it is the sum of conservative func- 
tional of the hydrodynamic energy J (p u2+ ™ 2 + pgz) dCl, of the mass JpdQ, and of the con- 

Q Q 

r r 2 2 

servative functionals J p lnp dD, and J q(z) dQ. Nonnegativity of the term p " ^ is obvious at 

n o 

p > 0. To prove nonnegativity of remaining summands of the integrand, we consider 

a = pgz+gH pin — 7— )+ (po(z)-p) 



We substitute p = po(z) (1 + if}- To be sure, the requirement p > is equivalent to rj > -1. For 
derivative with respect to variable rj, it is valid: 

^=Hln(l + rj)p (z)g, |^ = H^ypo( Z )s >0 at >j > -1. 

Unique minimum of the function a(n) corresponds to tj = 0. It is enough to be convinced that in 
the limit of small-amplitude waves the functional under study turns into ( 3Jj . □ 
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We term the functional (3.3| as a generalized functional of wave energy because in a small- 



amplitude limit it turns into wave energy functional (3.1 1 for linearized equations. At lack of 



wave disturbances, the integrand in 3.3| is equal zero. The idea of functional (3.3| form, is taken 



from the Hamilton formalism explained in |4|. Nevertheless, nonnegativity of the integrand, also 
some other important mathematical properties of the functional are discovered anew. 



Let us generalize (3.3| to a case when a solution may be nondifferentiable. The function- 



al of mass and hydrodynamic energy from which the functional is constructed (3.3|, should be 
conserved for nondifferentiable solutions; it follows from their physical sense. However the func- 
tional 

\ \ §H (p hx( — To))) is not conserved in a case if the density field is nondifferentiable. There- 
fore, we demand that the following inequality have to be fulfilled 

dH non i 



dt 



<0. 



(3.4) 



The sign "<" is chosen because the sign ">" in (3.4 1 leads to an unstable problem. We will use 



(3.3| and (3.4 1 as a basis of the theory of generalized solutions to nonlinear equations. Inapplica- 
bility of Hamilton formalism for considering of solutions with discontinuities of density follows 
from (33), pT4l >. 



3.1.3 Generalized solutions to the nonlinear problem 



Preliminary analysis. The equations (2.1 1 are written in a conservative form, in the form of 



fundamental conservation laws. From the physical point of view, the set of equations (2.1) is 
incomplete because it does not include the energy conservation law. Let co C Q be an arbitrary 
star domain with piecewise smooth boundary, and let S(a>) be a boundary surface of the domain 
a), dS = dSn, where dS is a surface element of S(a>), and n is an outer normal vector to S(co). We 
write an integral relation for energy 



J~e(x,z, t)dQ - ^ e{x,z, 0)dD. - J~ (|) e{x,z, t)v(x,z, t)dSdt = 0, e = p| 

w w o S(tu) 



7 7 
U + W 



+ g Z \ 



(3.5) 



If the solution of our problem is differentiable, then ( 3.5| follows from (2.1 1; but if the solution 
does not belong to a class of differentiable functions, then relation ( 3.5| does not follow (2.1 
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Therefore, permitting nondifferentiable solutions, we have to include (3.5) into the system of 
general equations as an individual equation. 



Eqations (2.1), (3.5| are not enough for statement of a correct problem. For stability of the 



nonlinear problem, relation (3.4 1 should be fulfilled as well. As conservation of energy functional 
( 3.5| is postulated, the inequality 

d 



dt 



(/(p)) + V(/(p)v) = v<0, /(p) = pln(p), 



(3.6) 



follows ( 3.4[ ) and ( 3.5| . Apparent requirement of compatibility of the relation ( 3.6| with the con 



tinuity equation (2.1 1 results in conclusion that v may be nonzero only at the points where p is 
nondifferentiable. 

Further we shall use / (p) = pln^^^y), where po(0) is the maximum value of the density (on 
the bottom), 

use of the function / of fixed sign is more convenient. 



Definition of a generalized solution. Using the set of equations (2.1 ), we routinely build some 
integral relations for definition of a generalized solution 



J~p s (x,z, 0)p(x,z,0)dQ + J" |p~J^ + P u ~^ + P w ~^f^ = ®' 

Q Q T 

J" u(x,z,0)p(x,z,Q)u s (x,z,0)dCl + ^ w{x,z, 0)p{x,z,0)w s {x,z, 0) dQ.+ 
n o 



(3.7) 



du s (x,z,t) 2 du s (x,z,t) 



du s (x,z, t) 



+ | \pw 

Qt 



-r h PUW ^— 

ox dz 

Qt ' 

dwAx,z,t) dwAx,z,t) ? dwAx,z,t) 

' +puw — S V +pw 2 — s \ ' +pgw s (x,z,t) 
dx -'- 



dt 



dz 



dQ T = 0. 



The relations should be fulfilled for any p s (x,z,t) e C^°(Q) x C°°[0, T], u s (x,z,t) = w s {x,z,t) 



~^t> 4>s(*>z>t) e C^ o (Q)xC oo [0,T], p s {x,z,T) = ip s {x,z,T) = 0. In equations (3.71, Q T = Q x [0,T]. 
For the basic functions (velocity components and density), it is convenient to introduce the stream 
function, u(x,z,t) = Jti*0l ) -w(x,z,t) = - ^'^' t \ where ip is a curve in W^iO.) parameterized with 
time t [25]], MM, EH- 
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However, the above integral relations are not sufficient for unambiguous definition of a gen- 



eralized solution. To be sure, we should supply the equations (3.7 1 with the relation (3.5) and the 



equation (3.6) in an integral form: 



^ f(p(x,z,t)) dQ- J~ / (p(x,z,0)) rfQ+J^ (j) f (p(x,z,t))v(x,z,t) dSdt = ^ ^ v (x,z, t)dDdt < 0. 

OJ co S(co) co 

(3.8) 

Here S(co) is a surface of star domain co with smooth boundary. Below, we formulate working 
definition of the generalized solution. 

We construct here a definition of a generalized solution, which should be acceptable from 
physical point of view and mathematically correct. From the previous reviewing, it is evident 
that ( 3.7[ ) is not enough to find uniquely a physically justified generalized solution. It is necessary 
to supply ( 3.7| with requirements ( 3.5| , ( 3.8| , or with equivalent ones. If 



p(x,z,t), p(x,z,t)u 2 (x,z,t), p(x,z,t)w 2 (x,z,t) g Lj (Qt)> 
p(x,z,0), p(x,z,0)u 2 (x,z,0), p{x,z, 0)w 2 (x,z, 0) e L 1 (Q), 



(3.9) 



then relations (3.7 1 have clear mathematical sense. Requirements (3.9) seem to be natural. There 



fore, at analysis of (3.5), (3.8), we assume that requirements (3.9) are fulfilled. If the solution is 



piecewise-continuous and restricted in Q T , then all integrals in (3.8), (3.5) exist; and it follows 



(3.81, that (3.51 



Je(x, z, t) dD. = e{x,z,0) dCl. 
n Jn 

^ f(x,z, ti)dCl < J~/(x, z, t 2 )dD. for all ti>t 2 , 
n n 



(3.10) 
(3.11) 



A hypothetical case when requirements (3.9) are fulfilled but the solution belongs to a class of 



functions, for which the surface integral in (3.5) is not defined, deserves special examination. In 



this case, it is possible to try to interpret a surface integral in (3.5) as a pair of functionals ( x and 
z components). It is important that these functionals exist for any star-shaped domain co with 
smooth boundary S(co) and the functionals are additive ones. That is, if co = coi U co 2 , oo\ n co 2 = 0, 
then a functional over surface S(co) is equal to sum of functionals over surfaces S(coi) and S(co 2 ), 
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and relations (3.10) and (3.5) are fulfilled. One can easily be convinced that if requirements (3.9) 



and (3.10) are fulfilled, then such functionals can always be defined. Therefore, for definition of 
the generalized solution, we demand realization of (3.9), (3.10); at the same time, relation (33) 



is not utilized in explicit form. Analogous analysis reveals that condition (3.11 1 has to play an 
important role in definition of the generalized solution, but it is possible to not use ( 3.8| resulted 



in (3.11). 



p{x,z,t) 



Definition 1. Let v(x,z,t) = (u(x,z,t),w(x,z,t)). We term X(x,z,t) ■ 

solution of equations (2.1), if u(x,z,t) = , u>(x,z,t) = - g^' > ^l^n = and z/(|3.7|) (|3.9|), 



v{x,z,t) 
I do 



as a generalized 



(3.10), (3.11), are fulfilled for any p s (x,z,t) e (Q) x C°°[0, T], u s (x,z,t) = h^, w s (x,z,t) = —$± 



ip s (x,z,t) G C M fQ)xC M [0J], p s (x,z,T) = ip s (x,z,T) = 0. 

Remark 1. For Navier-Stokes equations, the energy conservation law is not included into the system of 
constitutive equations, but the solution belongs the class V^lQr) 122]/; and the energy relation follows 
Navier-Stokes equations for functions belonging to V2(Qt) /132|/. j\33\]. In our case, an energy relation 



does not follow from integral relations ( |3.7 ), and the energy relation is imposed independently. Per se, 
we apply the condition of conservativeness (conservation of mass, energy) widely used in gas dynamics. 
In fact only the condition ( |3.I 1 ) important for inhomogeneous liquid, is novel. The condition ( [3.1 1| ) is 
important for viscous liquid as well. Fluid of constant density p = 1 was studied in known monographs 

EM, EM- 



4 Simulation of vortex destruction 



We take the initial condition: 



p(x,z,0) = p {z), ip(x,z,0) = A expi - 



,2 / \2~{ 
X-Xq\ [Z-Zq 



L 



(4.1) 



where l x = 0.16 m, Z z = 0.052 m, A = -0.0095 m 2 /s, x = 0.5 m, z = 0.125 m. This initial condition 
approximately corresponds to the disturbance created from an oscillating paddle in the tank in 
MackEwan's experiments [9j], [I10IJ . The horizontal scale is approximately twice more than the 
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vertical one. The amplitude of a vertical velocity is of 5 cm/s, the amplitude of a level velocity is 
of 16 cm/s. The top speed of internal wave propagation in the tank within a linear approximation 
is c = 10 cm/s. That is, the fluid velocity amplitude exceeds the limit of linear wave propagation 
speed. 

The solution procedure was grounded on numerical integration of the system of Euler equa- 
tions ( 2.1[ ) by finite-difference method. The fluid density for t = 3 s, 7 s, 8s is shown in fig. [T][2}[3j 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 



Figure 1: Fluid density at t = 3 s. H = 6.23 m, h = 0.25 cm. 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 



Figure 2: Fluid density at t = 7 s. H = 6.23 m, h = 0.25 cm. 

The overturning and formation of small-scale structures is advancing by formation of a tongue 
of a heavy liquid, penetrating into strata of a light liquid. Some layer structure containing seg- 
ments with inverse density is generated as a result. At t = 5 s, abruptions arise in tongues, and 
isolated fragments of a light liquid arise inside a heavy liquid. At t = 7 s, the process of breaking 
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0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 

Figure 3: Fluid density at t = 8s. H = 6.23m, h = 0.25cm 

becomes intensive. At t = 9 s, about 10% of fluid in the tank is retracted in intensive small-scale 
convection. The vortex is atomized and the square of the convective spot slowly increases. How- 
ever, with the course of time, small-sized blobs of a heavy liquid subside downwards, and blobs 
of a light fluid go upward. As a result, stable stratification is practically restored to the moment 
t = 14s. 

As a whole, the picture of wave breaking qualitatively coincides with the one circumscribed 
in [9]], [I10IJ . However, the last stage of reestablishing of continuous stratification is absent. It is 
explained by ignoring dissipative effects in the model. 

The flow function is shown in figs. |4}[5] At the first stage, the disturbance behavior is typical 




Figure 4: Lines of the flow function at t = 3 s. H = 6.23 m, h = 0.25 cm 

for internal waves and is perfectly explained by the theory of small-amplitude internal waves: 
left-hand and right-hand waves arise from this vortex. Because of the initial density perturbation 

14 



980 
976 
972 
968 
964 



L 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



Figure 5: Lines of the flow function at t = 7 s. H = 6.23 m, h = 0.25 cm 

is equal zero, the field of density perturbation is antisymmetric, and the field of a flow function 
is centrally symmetric. The symmetry is maintained in good approximation even when irregular 
movement is developed. 

Irregular structures in the flow function appear later than in the density. Perhaps, it is ex- 
plained by the fact that a flow function is an integral of velocity, hence it is the most smooth of 
considered physical fields. Two opposite jet flows are formed at t = 3 s, and then they in unison 
create considerable velocity shift. The energy of small-scale waves is scooped from a consider- 
able kinetic energy of the vortex due to instability development. It explains great intensity of the 
developed small-scale convection. 

To verify exactness of our simulation, the outcomes of simulation with larger spatial step 
h = 0.5 cm are shown for comparison in fig. [6j [7j 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 □ 960 

Figure 6: Fluid density at t = 3 s. H = 6.23 m, h = 0.5 cm 
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Figure 7: Fluid density at t = 7 s. H = 6.23 m, h = 0.5 cm 

At t < 3 s, the simulation with the grid step h = 0.5 cm coincide with the simulation with the 
scale h = 0.25 cm. At t = 5 s, some small-scale structures arise, and simulations with different 
steps differ a little. The distinctions increase with time. However, rather high concurrence of the 
outcomes is kept. All large-scale details coincide. 



4.1 Study of dependence of a vortex destruction on a stratification scale 

Identical simulations of a vortex destruction for stratifications with doubled and twice dimin- 
ished H have been carried out for study of dependence of the phenomena on a stratification 
scale. Some outcomes for H = 3.1 m are shown in fig. [8j [9j We see that when we decrease H, the 
evolution becomes an usual wave process. Fluid movement becomes more horizontal. The effect 
of wave breaking is starting later in spite of the fact that the inverse Brunt-Vaisala frequency is 
less. The wave breaking develops more slowly, and destruction runs languidly. Regular wave mo- 
tion of the fluid is maintained as a whole, and only separate "winged nuts" are generated. In fig. 



10] [TT] some outcomes of simulation of the same wave, but for stratification with doubled value 
of H, are shown. Comparison between the simulations relating to normal value of H and the 
simulations relating to diminished twice H displays that the phenomenon of vortex destruction 
and generation of small-scale convection is increased as H. 

Therefore, we see that the vortex destruction effect grows with H, and a case of almost homo- 
geneous fluid may be of interest. To consider the almost homogeneous fluid, we have carried out 
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Figure 8: Fluid density in the tank at t = 3 s. H = 3.1 m, h = 0.25 cm 




Figure 9: Fluid density in the tank at t = 7 s. H = 3.1 m, h = 0.25 cm 




Figure 10: Field density in the tank at t = 3 s. H = 12.4m, h = 0.25 cm 



17 




0.05- 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 



1000 

998 

996 

994 

992 

990 

988 

986 

984 

982 

980 

978 

976 



Figure 11: Fluid density in the tank at t = 7s. H = 12.4m, h = 0.25cm 

test simulation of evolution of the starting vortex for a stratification with very large H = 31 1, 5 m. 
The considered fluid is almost homogeneous and the density is varied only by 0.08 %. The simu- 



lation outcomes are shown in Fig. 12 13 14 We see the wave collapses, and the time of starting 



of wave breakdown is the same as in case H = 12.4 m. Therefore, for large H, the time of wave 
breakdown is determined not so much by value of H, as by starting conditions. If we consider the 
flow function, we discover the effect of formation of oppositely directed jet flows with velocity 
shear between them. Due to appearing this velocity shear, the instability arises. Small-scale fluc- 
tuations appear in the density field, and with some retardation, they appear in the flow function. 
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Figure 12: Fluid density in the tank at t = 3 s. H = 311.5 m. h = 0.25 cm. 



Some late instant of evolution of the flow function for the same initial vortex, but in the 
fluid of a constant density, is shown in fig. 15 The fluid density is constant and is equal to the 

18 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 



1000.0 
I 999.9 
999.8 
999.7 
999.6 
999.5 
999.4 
999.3 
999.2 
999.1 
999.0 
998.9 



Figure 13: Fluid density in the tank at t = 6 s. H = 31 1.5 m. h = 0.25 cm. 




typical water density. We see that the vortex breakdown and formation of small-scale convection 
is absent in a strictly homogeneous medium, and only the vortex has become more symmetric. 
Perhaps, it is explained by the fact that any function ijj(R), R = ^f(x 2 + z 2 ) gives stationary solution 
to equations for homogeneous fluid. Therefore, the vortex takes a symmetric form in time. 

Numerical simulations reveal that the phenomenon of breakdown becomes more and more 
brightly with increase of stratification scale H. Nevertheless, the phenomenon of vortex destruc- 
tion is absent, when the fluid density is strictly a constant. It points out that continuous limit 
from a stratified fluid to the case of a fluid of constant density is absent for our initial conditions. 
The approximation of constant density is out of the physical sense in such cases, even if density 
variations are very small. 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



Figure 15: The flow function lines at t = 7 s in a homogeneous fluid. 



5 Explanation of the effect 

Below we try to explain why evolution of a vortex in fluid of strict-constant density qualitatively 
differs from evolution of a vortex in fluid whose density is almost constant and varies within very 
narrow range. 

Let us take p = 1 for simplicity. Then the set of equations follows from ( 2.1| 



-=- + u^— + w^— = 0, (5.1) 

dt dx dz v ' 
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In (|5.1|), 5 = |f - fn> u{x,z,t) = w{x,z,t) = ~^r-. We assume that £(x,z,0) is a continuous 



<9t/> 



function, and max(|£(x,z,0)|) < M. At the first step, we formally write a solution to (5.1 1 using a 



characteristic method and supposing the functions u (x, z, t), w{x, z, t) being known. The character- 
istic method gives us that \£{x,z, t)\ < sup(\£(x,z, 0)|) and £(x,z, t) e L p (Q) for any p > 0. Applying 
the corollary of Sobolev embedding theorem [130] , [131] , we arrive at conclusion that the function 
ip(x,z,t) is continuous with respect to spatial variables, and u(x,z,t), w(x,z,t) are bounded ones. 
Moreover, the characteristic method also shows that all evolution of the starting vortex is only 
vortex deformation, but vortex destruction is impossible. 

In case of weakly stratified fluid, the equation for E, is like this 



d£ d£ 
dt dx dz 



(5.3) 



pQ = g 



dp dp 
dx dz 



du du du 
dt dx dz 



+ 



dp 

dx 



dw dw dw 

-zr- + U h W^— 

dt dx dz 



The equation for xjj is the same, as ( 5.2[ ). The right-hand terms containing square brackets are 
small against the first right-hand term in Q. Therefore, we ignore them and we consider only 
the approximation pQ « We discussed at beginning of the paper that a starting vortex in 
a stratified fluid by all means shapes tongues of light and heavy fluids which become very thin 
in time. Therefore, very thin density structures arise in time. So, after some evolution, the term 
\g^\ becomes essential. This source Q is of small scales and the one generates new vortices 
of smaller scales. These new small-scale vortices arise in weakly stratified fluid, and they should 
generate vortices of smaller scales in due time by means of the mechanism they arisen themselves. 
Therefore, we have a cascade process of generation of vortices of more and more small scales. 



Example 1. Any function ip(R), where R = Vx 2 + z 2 and ip{R) vanishes at infinity, satisfies ( 5.l|5.2 ). 
It gives a stationary solution. For example, we can take 



ip(R) = -M (K + L)exp(-*), 



(5.4) 



where L is a positive constant. We can easily do some analytical evaluations this such a function ip. This 
current function describes rotation of the fluid with velocity v(R) = = -yf f exp (-R/L) dependent 
on the distance from the vortex center. This velocity field has maximum at R = L. We can consider ( |5.4| ) 
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as approximate solution to (5.3), when we neglect of the term pQ ~ g^f-. Then we substitute (5.4) in 



continuity equation and calculate 



p (x, z, t) = p z - R * sin 



'v(R)t' 
2^R 



/ / fe exp(-R/L)t 



(5.5) 



At the following step, we calculate the source 



n d P ,1 P • / [gexp(-R/L)t 
pQ~g2 x =gp [ Z -R*smU r 



x . I gexp(-R/L)t\ I gexp{-R/L)t 

-— *sm \ I r — + cos /- — 

R VI 2nL VL 2nL 



2nL 
g exp {-R/L) tx 



(5.6) 



L 2nL 

Let stratification be linear and ^ = -a. We see that spacial derivatives of p increase in time up to 
infinite values, and it is the reason we have used a technique of generalized solutions. For small t we 



can neglect of the terms Uj^ + w^ in ( 5.3) and integrate the equation ( 



5.3) 



R\ i R" 



agx 
~RL 



2n(R + L) 



cos 



2tt 



-1 



exp 



iKM-f) 
f«*p(-f)'i 



+ i^^sin 



2n 



(5.7) 



Here corresponds to (5.4), and E,\ describes generation of new vortices. The expression for £i contains 



growing term - ^j^- sin | SXP 2 7t r j> an d due to presence of sin | j this term starting some t is a 

function changing sign as a function of coordinates. It displays generation of new vortices in time owing 
to stratification influence. 



We also note that other right-hand terms in Q in (53) take into account inertial forces arising 
in flow moving with acceleration. This inertia forces are equivalent to presence of some gravity 
field. That is, the effect of cascade vortex destruction in a nonuniform fluid can start without 
gravity field, but when we apply accelerated movement to vortex motion. 



6 Basic outcomes and conclusion 



The system of Euler equations for an incompressible stratified fluid was studied. A non-negative 
nonincreasing functional extending wave energy functional of the theory of linearized Euler 
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equations is suggested. The functional value is conserved over differentiable solutions, and di- 
minishes if density discontinuities are present. Functional properties have allowed using it for 
analysis of statement of a generalized problem. It is shown that Euler equations are not enough 
for statement of a correct generalized problem in case of a stratified fluid in the gravity field. 
Some auxiliary conditions are formulated and justified. One of them is an energy conservation 



law, and the second is a special requirement for density (3.6 1. The statement of a generalized 
problem is formulated. The finite-difference problem-solving procedure is developed. Existence 
of a weak solution is proved. 

The phenomenon of destruction of starting vortex in conditions of McEwan's laboratory ex- 
periments [9|> [HOI] has been simulated and studied. Qualitative concurrence of the simulated 
effect with the one observed in laboratory experiments, is good, but the last stage of restoring 
of smooth stratification is absent due to usage of an ideal liquid model. Some new data about 
the phenomenon, difficult for laboratory measurement, are obtained, and evolution of the flow 
function is studied. It is shown that large gradients appear in the velocity field at the expense of 
nonlinear effects and they may play an important role in development of instability. 

By means of numerical experiments, dependence of the solution on a stratification scale H 
is studied. It is revealed that the effects of vortex destruction and formation of small-scale con- 
vection grows with H for strongly nonlinear waves. On the contrary, when we decrease H, the 
effect of wave destruction becomes weak. The wave breakdown is starting from the density field 
and some later formation of small-scale structures in the flow function is observed. The effect of 
vortex destruction is not observed when the fluid density is strictly a constant. 
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